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Abstract 

Energy-momentum and angular momentum carried by electromagnetic field of 
two point-like charged particles arbitrarily moving in flat spacetime are presented. 
Apart from usual contributions to the Noether quantities produced separately by 
particles 1 and 2, the conservation laws contain also joint contribution due to the 
fields of both particles. The mixed part of Maxwell energy- momentum density is 
decomposed into bound and radiative components which are separately conserved 
off the world lines of particles. The former describes the deformation of electromag- 
netic clouds of "bare" charges due to mutual interaction while the latter defines the 
radiation which escapes to infinity. The bound terms contribute to particles' indi- 
vidual 4-momenta while the radiative ones exert the radiation reaction. Analysis of 
energy-momentum and angular momentum balance equations results the Lorentz- 
Dirac equation as an equation of motion for a pointed charge under the influence of 
its own electromagnetic field as well as field produced by another charge. 

PACS numbers: 03.50.De, ll.lO.Gh, ll.30.Cp 



1 Introduction 

The most natural and widely accepted equation of motion for a charge when radiation 
reaction is taken into account is the Lorentz-Dirac equation [1]. This equation has been 
discussed mainly for the case of one charge in an external electromagnetic field. In the 
present paper we consider an isolated system of two point electric charges and their 
electromagnetic field. We study the electromagnetic energy-momentum and angular mo- 
mentum radiated by charges; the study of energy-momentum and angular momentum 
balance equations implies the Lorentz-Dirac equation for more than one charge. 
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The dynamics of the entire system is governed by the action 



^ = E (-^« / + / dSaAa,^zli^ -^j d^vU^F", (l-l) 

where /^^ = '^ai^n^a,!^ — d^Aa,n) is the total field generated by two charges. Charge 
moves on a world line £ M4 described by functions z^{sa) which give the particle's 
coordinates as functions of proper time «„; :— dz^/dsa is the a-th four- velocity. 

The action is invariant under the space-time translations and rotations which con- 
stitute the Poincare group. This immediately implies conserved quantities which place 
stringent requirements on the dynamics of the system. They demand that the change in 
electromagnetic field momentum [2], 



p:^ = ]^da^T'^\ (1.2) 

and angular momentum [2], 

MZ = I da« {y^T^^ - y^T^^) , (1.3) 

should be balanced by a corresponding change in the total momentum and total angular 
momentum of the particles. (By dua we denote the vectorial surface element on a space- 
hke hypersurface E.) 

Since the Maxwell energy-momentum tensor density, 

is quadratic in the field and this field satisfies the superposition principle, the total elec- 
tromagnetic field stress-energy tensor is 

T^'' = T^;:^+Ti;^ + TZ, (1.5) 

where a-th particle density T^'^'^ is given by the expression (1.4) where "total" field 
strengths f^'' are substituted by "individual" ones f^^^y The mixed term 

^^TZ = fi;,]fk^ + fti^fh>^ - (/ri)/i? + /(2)/i;^) (1-6) 

describes the joint contribution due to both fields. 

In this paper we study radiation produced by an isolated system of two point electric 
charges and their electromagnetic field. Outgoing electromagnetic waves remove energy, 
momentum, and angular momentum from the sources which then undergo radiation re- 
action. The verification of conservation laws is not a trivial matter, since the interference 
contribution (1.8) involves divergent terms. 

In the derivation of particle's equation of motion, Dirac [1] evaluated the fiux of 
electromagnetic energy-momentum over a narrow world tube surrounding the particle's 



2 



world line. The author substituted the components of the retarded Lienard-Wiechert 
field in the stress-energy tensor (1.4) inside the narrow tube. In 1970 Teitelboim [3] 
splits the "retarded" stress-energy tensor into the "bound" and "emitted" parts which 
are separately conserved off the world line of the particle. The author calculates the 
flow of energy-momentum out of the portion of Bhabha world tube [4] bounded by tilted 
spacelike hypersurfaces which are orthogonal to particle's four-velocity at instants r and 
T + dr, respectively. Bound part, Tbnd, describes a rigid electromagnetic "cloud" which 
are permanently attached to the source and carried along with it. "Bare" charge and 
"cloud" constitute new entity: dressed charged particle. Tbnd contributes into particle's 
inertia: 4-momentum of dressed charge contains, apart from usual velocity term, also a 
term which is proportional to the square of charge e^, 

P^ = ^a<-^<. (1.7) 

(Coulomb-like infinity stemming from the pointncss of "bare" source is absorbed by the 
rest mass within the renormalization procedure.) Time derivative of the second term 
in eq.(1.7) is the well-known Schott term which describes a reversible form of emission and 
absorption of field energy, which never gets far from the point-like source. The radiative 
part. Trad; yields the Larmor relativistic rate of radiated energy-momentum. It detaches 
itself from the charge and leads an independent existence. This rate together with the 
Schott term constitutes the Abraham radiation reaction vector. Lopez and Villarroel [5] 
split the torque of the stress-energy tensor into bound and emitted components which 
possess analogous properties. 

The results can be applied to the first and the second terms in the total electromagnetic 
field stress-energy tensor (1.5) which describe "individual" radiation contributions due to 
particles 1 and 2, respectively. The question is what part of the mixed density should be 
taken instead of (1.6) to describe the radiation which reaches to a very distant sphere? 

Aguirregabiria and Bel [6] studied the interference part of energy-momentum, 

PL = ld<^,TZ, (1.8) 

carried by electromagnetic field of two point charges. The authors prove the fundamental 
theorem that the mixed radiation rate does not depend on the shape of spacelike surface 
E which is used to integrate the mixed part (1.6) of the Maxwell energy-momentum 
tensor density. For a prescribed plane motion of the charges the perturbation scheme is 
elaborated within the framework of predictive relativistic mechanics [7, 8]. The lowest 
approximation gives the well-known expression [9, p. 2 14] for the dipole radiation of two 
point charges moving according to Coulomb's law. In Ref. [10] the scheme is applied to 
the angular momentum carried by electromagnetic field of two point charges. 

In the case of particles we would merely obtain as an obvious generalization of 
eq. (1.5) the sum of A^ one-particle terms and the mixed contributions corresponding to 
N[N—l)/2 pairs of charges. Hence, the interference component dominates in the radiation 
from a bunch of identical charged particles (e.g., in free electron lasers). To evaluate the 
radiation of a relativistic A^-body system, Klepikov [11] defines the center of a system of 



3 



radiation events which allows to synchronize the instants at which electromagnetic waves 
emitted by different charges combine on a very distant sphere. Fourier analysis is applied 
to calculate the time and angular distributions of energy-momentum flux. The radiation 
of a bunch of charged particles moving in a uniform magnetic field is considered in detail. 

The only exact solution is obtained by Rivera and Villarroel Ref. [12, eq.(3.27)]. 
The authors calculate the rate of radiation (including interference part) generated by two 
identical point charges rotating uniformly at opposite ends of a diameter, in a fixed circle. 
External fields which govern the strictly prescribed motions are constructed. The rate of 
radiation is evaluated via the retarded Lienard-Wiechert fields produced by the charges. 
Further [13] more general case of circular motion of two unlike charges in two coplanar 
and concentric circumferences is considered. 

Note that nine years before the radiation by a system of two uniformly circling charges 
has been evaluated by Hnizdo [14] (see also discussion [15, 16]). The author concludes 
that "the power radiated by such a system equals exactly the rate at which work is done 
on the system by external force" . 

In this paper we study the interference part of energy-momentum and angular momen- 
tum of electromagnetic field generated by two arbitrarily moving charges^. We restrict 
ourselves to the retarded Lienard-Wiechert solutions; the advanced ones are rejected on 
the grounds of causality. In Section 2 we introduce coordinate system which makes rela- 
tively easy the calculation of the covariant 4-momentum radiated by interacting charges. 
The calculation is performed in Sections 3 and 4. We reveal divergence-free radiative 
part of the mixed density (1.6). It determines the radiation that escapes to infinity while 
the (short-range) bound part modifies individual 4-momenta (1.7) of dressed particles. 
The mixed part of radiated energy-momentum depends only on velocities, accelerations 
and the relative 4-position of the charges. In Section 5 we derive equations of motion 
of interacting charged particles. Analysis of energy-momentum and angular momentum 
balance equations results the well-known Lorentz-Dirac equation. In Section 6 we study 
symmetry properties of radiative energy-momentum and angular momentum which rely 
on invariance of action (1.1) under inversions of space and time axes. In Section 7 we 
discuss the results and implications. 

2 "Interference" coordinate system 

To perform the surface integration (1.8) of interference stress-energy tensor, an appropri- 
ate coordinate system is necessary. Such a coordinate system is introduced in Ref. [6]. 
It involves the evolution parameter A associated with an inertial observer; the surface of 
integration is a surface of constant A. In Refs. [17] and [18] this coordinate system is 
adapted to the simplest hyperplane T,t = {y E M.4 : = t} associated with an unmoving 
inertial observer. The "laboratory" time i is a single common parameter defined along 
all the world fines of the system. 

^The generalization of this work to N charges is an obvious one. 
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Figure 1: The interference picture in a plane S^. Sphere Sa with radius = t — ta is 
centered at point Oa with coordinates {to), i' = 1, 2, 3. Angle (p distinguishes the points 
of intersection 81(182 — C{0, h) which constitute support of integrals (2.1) and (2.2). 
/c^, and h arc the components of the future oriented null 4-vector A;" = a(y° ~ ^ai^a))- 
Matrix Vt determines transition to "momentarily rotating" Lorentz frame. 

The mixed contributions to energy-momentum, 

P:^t(^) = / dao7;°^ (2.1) 

J St 

and angular momentum, 

MZit) = / dao {y'^Tf:, - y'^ltt) , (2-2) 

^St 

are due to interference of spherical wave fronts 81 and 5*2 in (see Fig. 1). Sphere, 

8a{Za{ta),t - Q = {?/ G M4 : (/ - ta)' = " 4(^a))', ^ t,t - ta > 0}, (2.3) 

i 

is the intersection of the future light cone with vertex at point Za{ta) G Ca and S^. This 
contribution is zero if the relative position 4-vector g = 2:1 — 2:2 is timelike. If q is spacclike, 
the intersection 5"! fl 5*2 is the circle C{0,h) with radius h; in "momentarily rotating" 
Lorentz frame the circle C{0,h) lies in Oxy plane and centered at the coordinate origin 
(see Fig. 1). If points zi and Z2 are related by a null ray, the intersection fl 5*2 contains 
the only point. 

2.1 Local map 

To find the local expressions for coordinate transformation (y") 1— > {t,ti,t2,<f), we trans- 
late the origin of the laboratory Lorentz frame at the center O of the circle C{0,h) — 
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Figure 2: The past light cone with vertex at point y E T^t is punctured by the world lines 
of the 1-st particle and the 2-nd particle at points Zi{ti) and Z2{t2): respectively. The 
vector Ka is a null vector pointing from the emission point Za(ta) = {ta, zl^ita)) to a field 
point y. The relative position 4- vector q — zi — Z2 is equal to difference K2 — Ki] its 
square {q ■ q) = — 2(i^2 • Ki). 



31(182 and then rotate space axes till a new z-axis be directed along 3- vector q := zi — Z2, 

y'^^z:{ta) + n\'{h,t2)k:'. (2.4) 
Here ka, a = 1,2 is the future oriented null- vector with components 

k'a-t-ta, kl^hsin^, kl^hcos^, kl^i-irl+^-^^^^^, (2.5) 

which arise from analysis of triangle O1O2H pictured in Fig. 1 (we denote q = |q|). 

Matrix space-time components are = ^^0 = ^tiO- Its space components Qij constitute 
an orthogonal 3x3 matrix (A. 6) which determines the rotation described above (see 
Appendix A). It defines new orthonormal basis, 

= cos (fq cos 'j?qei -|- sin (fg cos '&qe2 — sin i?qe3 , 
= - sin (fiqei + cos (pqe2 , (2.6) 
Uq — COS (pq sin i?gei -|- sin (pq sin i?qe2 -|- cos ■J^qBs , 

which is constructed from components of the relative position 3- vector q, e.g. cosf^g = 
?VvW + W, cos^, = gVq. 

To find the Jacobian of coordinate transformation (2.4), we derive its differential chart. 
Setting a = in eq.(2.4) immediately follows y° = t. Since t — y°, then dt/dy°' — Soa- 
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Because y and Zaitg) lie on the light cone (see Fig. 2), a change field point y comes with 
a change in ta- Suppose that y is displaced to the new point y + Sy. The new intersection 
of the past light cone of this vertex with the a-th world line is then Za{ta + Sta)- These 
points are still related by the equation 

(y° + 5y' -ta- 5taf = ^(2/^ + 5y' - + 5ta)f. (2.7) 

i 

Expanding this to the first order in 5y and 5ta and using the cone equation (2.3), we 
obtain Ka^aSy"' - {va ■ Ka)8ta = or 



Here Ka is a-th null vector pictured in figure 2 and symbol denotes the scalar product 
[va ■ Ka), taken with opposite sign; noncovariant 4- velocity Va := (1, dzydta). 
For the angular variable we have 



dy' 

where 



= n^'K ' (2-9) 



45 = 0. *i = ^, kl^-'^. fcj = 0. (2.10) 

Recall that h is the radius of the circle C{0, h) = SiD S2 pictured in Fig. 1. 

Determinant of the matrix which defines this differential chart gives the inverse Jaco- 
bian: =q/(rir2). The "interference" surface element, 

dao = ^dtidt2dip, (2.11) 
is ill defined if and only if the particles are very close to each other. 



2.2 Global mapping 

Setting a = 1 in eq.(2.4) we obtain the coordinate system centered on an accelerated 
world line of the first particle. The fiat spacetime M4 is a disjoint union of hyperplanes 
= {y ^ : y^ = t}. An interference hyperplane is a disjoint union of retarded 
spheres Si{Oi,t — ti) centered at points Oi e with coordinates Zi(ti). A sphere is 
covered by its intersections with spherical wave fronts of the second source. Each circle 
Si n S2 can be labelled by the individual time t2 of the second particle and each point on 
a given circle can be labelled by its polar angle (/?. 

Going along the world line Ci we arrive unavoidably at the point t\^^{t), such that the 
future light cone of Zi[fi^{t)\ touches the 2-nd world line at point Z2{t) G (see Fig. 4). 
Light cone of upper vertices do not intersect the second world line at all. 

In context with the principle of retarded causality, is divided into two regions where 
outgoing waves sourced by charged particles combine in quite different manner. 
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Figure 3: For a given ti the wave front y^i is covered by circles Si fl S2 if the parameter 

t2 increases from tf^{ti) to t2'^^{ti). Minimal value labels the vertex of forward light cone 
which is punctured by (i at Zi(ti). The largest sphere touches the sphere Si at point 
N. World line ^2 punctures the future light cone of zi{ti) at point 2:2(^1'^^). Intersection 
of 5"! and the smallest sphere S^^ contains the only point S. 



(i) Causal, which is filled up by spheres 5"! of radn larger than or equal to t — t[^^(t). 

(ii) Acausal, where parameter ti increases from t\^^{t) to the instant of observation t. It 
is the ball bounded by sphere of radius t — t'i"\t) centered at point zi[t^^*(t)]. 



2.2.1 Causal region 

Causal region is spanned by curvilinear coordinates (2.4) where ti increases from —00 
to the instant t'^i*{t). To cover the sphere Si where ti is fixed we change the parameter 
t2 which labels point 2:2(^2) € C2- The starting point is the solution t'^2*{ti) of algebraic 
equation g° = q or 

ti-t2 = q(ti,t2), (2.12) 

where points Zi{ti) G Ci and Z2{t^2*) G C2 are linked by a null ray. The largest sphere 
5*2(02'^*, t — ^2^*) touches a given sphere 5*1(01, t — ti) at only point N (see Fig. 3). If 
parameter t2 increases to ^^'"{ti) being the solution of algebraic equation g° = — q or 

i2-ii = q(ii,i2), (2.13) 

the intersection fl 5'2'^" contains the only point S. If parameter ^2 changes from i2^*(ti) 
to t'^'"{ti), the sphere ^i is covered by circles C(0, h) — SiH S2. 
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Figure 4: The instant is the solution of algebraic equation (2.13) with t2 = t. 

Sphere Si{z'[^^, t — tl^*{t)) bounds acausal region. For a given ti e [t'[^^{t),t] the sphere 5"! 
with radius t — is a disjoint union of circles 81(182 if parameter t2 increases from t2^*(ti) 
to t2{t,ti). The maximal value satisfies algebraic equation ^2 + ^1 = |q|. The sphere 5*2' 
with radius t — t'2 touches Si at point S. Intersection 5*1 Pi 5*2^* contains the only point N. 

2.2.2 Acausal region 

Acausal region of an interference hyperplane corresponds to the fragments of the world 
lines which are not related to each other. (By this wc mean that the radiation emitted by 
the first particle during the interval [t{'^^ (t) , t] does not come to the second one and vice 
versa.) Nevertheless, the outgoing waves of these portions of world lines combine in . 

Acausal region is filled up by spheres 8i{Oi,t — h) e E^, where ti e [t'['^^{t),t]. The 
sphere 81 with fixed ti is the disjoint union of circles C{0, h) = SiH 82 if the parameter 
t2 increases from t'^2*{ti) to ^2(^1^1)- The starting point of this interval is still the solution 
of eq.(2.12) while the maximal value of ^2 satisfies the algebraic equation A;^ + A;° = q or 

2t-ti-t2^qih,t2). (2.14) 

A given sphere 6*1(01, t — ti) touches 5*2(02) ^ ~ ^2) o^^Y Point 8 (see Fig. 4). 

2.2.3 Surface integration 

In an analogous way wc construct the coordinate system centered on the world line of 
the second particle. If ^2 G] - oo,t^"*(t)] then h E [t[^\t2) , tf" {12)]; if t2 G [tr*W>^] then 
ti e [tT^{t2),t'i{t, ^2)], ^ G [0, 27r[. The ends of intervals are defined implicitly by algebraic 
equations (2.12), (2.13), and (2.14). 

The surface integration (2.1) and (2.2) can be performed via the coordinate system 



9 



centered on a world line either of the first particle, 



I I ^^^^ / / / ' 



or of the second particle, 



dip 



q 



(2.15) 



. . t t[{t,t2) 

J dt2 J dti+ J dt2 J dti j > 



dip 



rir2 



(2.16) 



To calculate the flows (2.1) of the mixed electromagnetic fleld energy and momentum 
which flow across the hyperplane E^, we should integrate the Maxwell energy-momentum 
tensor density (1.6) over angular variable ip and over time variables ti and t2- 



3 Angular integration of energy- moment una and an- 
gular momentum tensor densities 

In this Section we trace a series of stages in integration of the mixed Maxwell energy- 
momentum tensor density over ip. In Appendix A we derive some useful expressions. 

In terms of Minkowski coordinates {y") the electromagnetic field generated by a-th 
particle is given by 

/(a) = -^Ua Aka + — [tta AK + {K ' aa)Ua A k„] , (3.1) 
•"a 

where symbol A denotes the wedge product. We use sans-serif symbols for the retarded 
distance^, 

ra^-Vap{y''-Z:{Sa))u^M, (3-2) 

and for the null vector Ka — y — Za{sa) rescaled by a factor r~^, 

K-^{y"-4{sa)). (3.3) 

'a 

To express field strengths in terms of curvilinear coordinates {t,ti,t2,ip), it is advantageous 
to replace the retarded proper time Sa{y) by evolution parameter ta- The components of 
particles' 4- velocities Ua and 4-accelerations Ua, a — 1,2, become [2] 

K = laVi^iQ, = ^t(Va ■ VaK + iX, (3-4) 

where 4- vectors = (1,^^(4)), vj^ = {0:Vl{ta)) and factor 7^ [1— v^]~-^/^. Substituting 
these into eq.(3.1) and using the relation = K^/^a yields 

ha) = e„ ^3 + j ' (3-5) 



^Because the speed of light is set to unity, r^, is equal to the spatial distance between 2a[s™*(y)] and 
y as measured in momentarily comoving Lorentz frame where u" = (1, 0, 0, 0). 
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where 

r„ = X°-(K,v„), c„ = 7-' + (KaVa). (3.6) 

Note that r<j is the retarded distance (3.2) rcscaled by a factor 7^, i.e. Va = Ja^^^a- The 
separation vector has the form Clka, where components of null vector ka are given by 
eqs.(2.5) and matrix Cl determines the transition to momentarily comoving Lorentz frame 

associated with basis (2.6). 

It is straightforward to substitute the components of electromagnetic fields (3.5) in 
terms of "interference" coordinates {t,ti,t2,(p) into integrands of expressions (2.1) and 
(2.2) to calculate the interference part of radiated energy-momentum and angular mo- 
mentum, respectively. Integration of the mixed stress-energy tensor over angular variable 
is the key to the problem. All (/7-dependent terms are concentrated in the following con- 
structions: 

f • = ^ /d*'— ■ = /d^^, (3.7) 

27r J qrir2 27r J qri(r2)^ 



2^ J ' ^ = ^ I 



q(ri)V2' 27: J q(ri)2(r2)^' 







They are labeled according to their dependence on the combination of components of the 
separation vectors Ki and K2: factor a is replaced by K1K2, K^, K^, or 1 for VI2, V^, 
T>2 or 1)°, respectively. (The others and A"' are marked analogously.) 

The mixed part of the stress-energy tensor (1.6) is symmetric in indices 1 and 2. 
Substituting (3.5) into the first term of this expression and using the identities K2—K1 — q 
and {K2 ■ Ki) = -l/2{q ■ q) yields 






^'""^ dt,dt2 ' ' dtldt2 ' ' dt,dt2 ' ' dt^dtl 

// „ da /, ,/ da _n / // „ da ,, „ da ,, ,, d'^a 

- - '^<"'^m; - ^ ("■ a^; + % + '"''^atlm., 

after addition of similar terms and integration over (/?. World function (7(^1,^2) of two 
spacelike related points, zi{ti) e Ci ^-nd Z2{t2) £ C2) is equal to one-half of the square of 
vector q — z\ — Z2.I taken with opposite sign, 

a{h,t2)^-ll2{q-q). (3.9) 

Each second order differential operator, 

= + + ^ V + (3-10) 

OT1OT2 dti Ot2 
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has been labeled according to dependence of coefficients (3.7) on the combination of 
vectors Ki and K2. 

For the convolution /("^/^J, we obtain 

27r 



^ / difJf^'^C = eie2r°(A), (3.11) 



where function 

d'^a da da 
A = (7^— — -— — 3.12 

OtiOt2 Oti Ot2 

depends on two-point function (3.9) and its derivatives in time variables. 

To distinguish the partial derivatives in time variables, we rewrite the operator (3.10) 
as the sum of the second-order differential operator, 

and the "tail". 

For a smooth function /(ti,t2) we have 

r«(/)=n«(/) + /7r". (3.15) 

Cumbersome calculations which are presented in Appendix A give the relations 

(3.16) 



7r° 


= 0, 






< 








''12 









+t'2(cr-^)-«^^ 

which allow us to rewrite the sum of integrals (3.8) and (3.11) in terms of differential 
operators 11" and partial derivatives in ti and t2, 

P!? = (/S/wc.-^/,°?/S') (3.17) 



2 t \dhdt2J ' V ^^1/ V 



dti \ ' 'dtidt2j dt2 V dtidt^, 
d_ 



12 



We denote V21 the integral over ip of the remaining terms involved in tensor (1.6). It can 
be obtained by interchanging of indices 1 and 2. 

Setting /J, — and u — i in eq.(3.17), we obtain the first term of the mixed space-time 
components of the stress-energy tensor (1.6). Wc add the term where indices 1 and 2 are 
interchanged. Since zeroth components and fc^ of the separation four-vectors Ki and 
K2 do not depend on ip, the final expression get simplified, 



int 



^12 "T ^21 



6162 

2 

d_ 



ni 

d_ 



(3.18) 



dti 



9^2 j^O 



where 



Ai 



, n 9(7 



, n 9(J 



(3.19) 



Similarly we derive zeroth component Setting fi = and u = in eq.(3.17), 

we obtain the first one-half of desired expression. The second one, V21, can be derived 
via interchanging indices 1 and 2. The integral of energy density over the angular 
variable has the form 



^int — ' 12 "T ' 21 



dtidt2 



9S 9S 
dfidh 



(3.20) 



where three-point function. 



E(i,ii,t2) = 2A;i'/c^ + (7(ii,i2), 



(3.21) 



depends on particles' positions referred to the moments ti and ^2 before observation instant 
t as well as on t itself. 

We now turn to the integration of the angular momentum tensor density (2.2) carried 
by the electromagnetic field due to two pointhke charges. We present the torque m^^^ — 
y^TP^^ - y^T^^^ in the following form: 



m 



int 



IIV I IIU VII VII 

^12 + ^21 - ^12 - ^21: 



(3.22) 



where 



m>i^2 = (^r + K^) ^ 



fOXfv ^ldvfal3A2) 
/(1)/(2)A J{l)Jal3 



(3.23) 



It is straightforward to substitute the fields (3.5) into this expression to calculate the first 
term of expression (3.22). The others can be obtained by interchanging of the pair of 
indices (1,2) and (//, i^). 
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Having integrated expression Jm^2 over ip we obtain 



6162 



12 



dXi 



dU 



1, <9'^Ai „„ 



dt2 ' 'dtidt2 
P2^ + r°(^^2'Ai)-^;KAiB° 



f-t (A) + r° (^'^A) - ynB"" - 

OT2 



(3.24) 



where functions A and Xa are given by eqs.(3.12) and (3.19), respectively. 

Usage of the equahties in eq.(3.16) derived in Appendix A allows us to rewrite the 
integrand (3.24) as follows: 



6162 



12 



+ m 



^^12 



dXi 



d 



,dXu 



dU 



(3.25) 



d_ 

d 



+ nV^r^Ai)-^ A: 



II <9Ai 
dt, 



d 



(A) + If (^f A) - ^ «AP°) 



The other terms of mixed angular momentum, 



12 



21 5 



(3.26) 



can be obtained via interchanging of indices (1,2) and (/x, v). 

Wc sec that the integration of the mixed stress-energy tensor (1.6) over (/9 yields the 
combinations of partial derivatives in time variables. In the next Section we classify them 
and reveal the long-range terms which contribute into radiated energy-momentum. 



4 Radiative parts of mixed energy- momentum and 
angular momentum 

In previous Section we integrate the mixed part of the stress-energy tensor and its torque 
over polar angle. Resulted expressions describe contributions to electromagnetic field's 
energy-momentum and angular momentum due to interference of spherical wave fronts 
of charges e\ and placed at fixed points z^it}) G Ci and z^i^'i) € C2, respectively. The 
crucial issue is that the integrals (3.18), (3.20) and (3.26) have the remarkable property 
of being the sum of partial derivatives in time variables. This circumstance allows us 
to calculate how much electromagnetic field's energy-momentum and angular momentum 
fiow across a hyperplane E^. 
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It is natural to integrate the expression being the time derivative with respect to t2 
according to the rule (2.15), 



dt2+ j dti j dh 
t 



dG2{h,t2) 
dto 



(4.1) 



= j dtiG2[U,tf''{h)]- j dhG2[h,f2'\u)]+ J dhG2[U,t'^{t,U)]. 

-oo -oo tl^^t) 

Having apphed the rule (2.16) to the expression of type dG^/dti, we obtain 



dt2 j dti+ j dt2 
t 



dti 



dGi{h,t2) 



(4.2) 



= j dt2G,[tf''{t2),t2]- J dt2G,[tY\t2),t2]+ J dt2G,[t[{t,t2),t2]. 

-oo -oo t2''*(*) 

The end points are valuable only in the integration procedure. They are solutions of 
algebraic equations (2.12), (2.13) and (2.14). The retarded instants fa^itf,) and advanced 
ones tl'^^{ta) label the points N and 5" in which fronts of outgoing electromagnetic waves 
produced by charges touch each other (see Fig. 3). All the moments are before the 
observation instant t, so that the retarded causahty is not violated. 

It is worth noting that the functions ti^*(t2) and tf'^''(ti) are inverted to each other 
as well as the pair of functions t'l'^'"{t2) and V^^iti). For a fixed laboratory time t the 
functions t[{t,t2) and t'2{t,ti) are inverses too. These circumstances allow us to change 
the variables ta ^ i'a^{h) in the "advanced" integrals in eqs.(4.1) and (4.2). Further we 
couple them with their "retarded" counterparts. Since 



dtr\t2) 
dU 



Vlllq 



dtl'\ti) 1 + (V2nq 



1 - (VsHq) ' 



dti H-(vinq)' 



(4.3) 



we obtain 



Q2 + Q1 



dti 



1 - (VlHq) 



Gi — G2 



-00 
t 



+ 



+ 



dt2 



1 - (V2nq) 

1 + (VsHq) 



z 

I 



dU 



1 + (VlHq) 



G, 



^ - lV2nqj 
1 + (VlHq) 



Go 



tl=tr*(t2) 



tl=t[{t,t2) 



(4.4) 
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The unit vector = q/q is the third vector of orthonormal triad (2.6). The last integral 
is due to interference of outgoing electromagnetic waves radiated out by the particles 
within the acausal region (see Fig. 4). We take into account that 

dt'2{t,ti) l + (vinq) 



dti 



1 - (vanq) 



(4.5) 



Of course, one can change the variables tb i— >■ t^'^^{ta) in the "retarded" integrals in 
eqs.(4.1) and (4.2) and add them to their "advanced" counterparts, 



—00 



+ 



dti 



1 - (v2n,J 

Lri — rtj2 

1 - (VlHq) 

G1 + G2 



-lti=tf^{t2) 



1 + (ViHq 



—00 
t 

I 



1 + (VsHq) 

G1 + G2 



1 + (ViHq 



1 - (V2nq) 



+ / dti 

A combination of the "retarded" and the "advanced" terms is valuable too. 
Integral of a mixed double derivative can be written in the form either 



dU 



dti 



ti(*'*2) 

dt2 / dti 



+ 



—00 
t 

I 



dU 



dto 



tr*(*2) 

dG(h,t2) 
dt2 

dG{hM) 
dU 



d_ 



dG{hM) 



tl=t[{t,t2) 



- / dt2 



<9G'(ti,t2) 
dto 



ti=tr*(*2) 



(4.6) 



(4.7) 



or 



J dti j dt2+ I dti J dt2 

-00 tY''{ti) 



d_ 
at^ 



dG{h,t2) 



dti 

—00 
t 



dG{h,h) 
dti 

dGjtiM) 
dti 



t2 (ii) 



t2=t'^{tM) 



dti 



dG{ti,t2 

ati 



t2=t5"*(tl) 



(4.8) 
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Figure 5: 



Function ^[ti, 



is associated with the interference of wave fronts 
Si{Oi,t - ti) and S2{0^^\t - t^"*) at point N e E^. Function G[ti,t'2{t,ti)] is con- 
nected with the combination of wave fronts Si{Oi,t — ti) and 5*2(02, t — ^2) point 
S G Tit- If ti ^ t, both the point and the point S tend to Zi{t) = (1(1 E^. If 
lim G[ti, t2^^{ti)] = hm G[ti,t2{t, ti)] the value of integral of mixed derivative d^G / dtidt2 

over time variables does not depend on the order of integration. 

The question is what expression should be used. 

To compare Ti and T2 we change the variables in the "advanced" integrals and subtract 
(4.8) and (4.7). We arrive at the integrals being functions of the end points only, 



T1-X2 = G[h,tl'\t,)] 



+ 



either G[t[{t,t2),t2] 
or - G[ti,t'2{t,ti)] 



G[tT\t2),t2] 



t2- 



(4.9) 



It vanishes if and only if (i) limiting values of G (that evaluated at the remote past) cancel 
each other, and (ii) function G is smooth at points at which the world lines puncture Ef. 
Indeed, the part of difference (4.9) which depends on the momentary state of particles' 
motion can be rewritten as follows: 



\im{G[hX2\t,)]-G[t,MtM)]] 



(4.10) 



]im{G[tl'\t2),t2]-G[t[{t,t2),t2]}- 

t2 — >t 



The situation is illustrated in Fig. 5. 
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4.1 Criteria 

The main task of the present paper is to decompose the interference part of the Maxwell 
energy-momentum tensor density into bound and radiative components. The former 
modifies individual 4-momenta (1.7) of dressed particles while the latter shows how a 
charge is influenced by radiation of another charge. 

To reveal meaningful radiative part of mixed energy-momentum (2.1) and angular 
momentum (2.2), we apply the criteria which were first formulated in Ref. [3, Table 1]. 

• The bound part diverges while the radiative one is finite. 

• The bound component depends on the momentary state of the particles' motion 
while the radiative one is accumulated with time. 

• The form of the bound terms heavily depends on choosing of an integration surface 
while the radiative terms are invariant. 

There are, however, a several properties that the desired expressions have possess before 
they can be accepted. We fist them below. 

1. Radiative parts should be completely determined by particles' motion; they can not 
depend on distance to a point of observation. 

2. They should be produced by divergence- free expressions. 

3. Non-accelerated charges do not radiate. 

4. Balance of Noether conserved quantities yields the Lorentz-Dirac equation. 

4.2 Radiative part of mixed momentum 

To decompose the momentum (3.18) into bound and radiative components is a straight- 
forward integration of all the terms over time variables. Scrupulous computations reveal 
candidate for radiative part. 



which produce the terms that satisfy Teitelboim's criteria. 

Equations (4.1) and (4.2) imply that its integral over time variables is completely 

determined by values of the arguments of time differential operators at the ends of in- 
tegration intervals. When a consideration is restricted to the end points where radius h 
of intersection 5*1 fl S2 vanishes, the coefficients (3.7) get simplified. Since eq. (2.5) the 




(4.11) 
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integrands do not depend on the polar angle at all. The integral over 93 becomes unit 
operator. Radiative momentum (4.11) contains also non-trivial constructions, 



(4.12) 



which are calculated in Appendix A. All the expressions (A. 35), (A. 41), (A. 43), and 
(A. 53) are regular if /i = 0. Values of relative position 4-vector g, retarded distances ra, 
and some basic functions at limiting points are collected in Table 1 (see Appendix B). 

The sum of differential operators and W2 contains mixed double derivative of the 
function. 



■v\. 



(4.13) 



dtidt2 UL1UL2 

A surprising feature of integration of expression (4.11) over time variables is that the 
result heavily depends on the order of differentiation in ff^G^ / dtidt2. If one choose the 
rule (4.7) they obtain 



pU2i= / dh-f^^G\,[tiA 



J —00 



dii7r'G'2i[^i,^2*(^i)]- 



(4.14) 



Having applied the rule (4.8) we arrive at 



^'rad,12 



dhl2'G\2[tT(t2)M] + 



dt2l2'G\2[tt\t2)M]- 



Function under integral signs. 



13 % ^ [qab-Vb) 



(4.15) 



(4.16) 



+ 



is referred to the retarded and the advanced instants. Relative distance 4-vector qat '■= 
Za — Zb] in our notation qi2 = q while ^21 = —Q- Denominator Vb = —{qab " Vb) is the 
retarded distance of type (3.6) where both the field point and the point of emission are 
placed on particles' world lines. So, in expression (4.14) the field point is zi{ti) e Ci 
where ei is placed. Points of emission are Z2[t2^^{ti)] £ C2 in the first path integral and 
^2(^2'^^ (^i)] G C2 in the second one. 



r^2'' = q [1 - (V2nq)] 



r«'^" = -q[l + (v2nq)] 



[ti,t5«*(ti)] 



(4.17) 



[ti,tf^{ti)] 



In expression (4.15) the field point is in 2:2(^2) G C2 where the second charge is placed. It 
is acted on by the first charge placed at points either Zi[tl^*{t2)] G Ci or Zi[t'f'" {t2)] G Ci- 
The distances are as follows: 



^ret 



q[l + (ViHq 



„adv 



-q[l - (ViHq) 



(4.18) 
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It is of great importance that the integrand (4.16) does not depend on the distances 
to points N or S on the observation hyperplane E^. It depends on the relative distance 
between charged particles, on their velocities, and on their accelerations. The situation 
looks like a-th charge is acted on by 6-th one directly. The charges are connected by a null 
ray: in the first integral in eq.(4.14) or in eq.(4.15) the interaction be forward while in 
the second one backward in time. All the moments are before the instant of observation 
t which labels the upper limits of path integrals. 

Since the function (4.13) is not smooth in neighborhoods of intersections Ca ^ E^, the 
expressions (4.14) and eq.(4.15) are not equivalent. By virtue of the relation (4.9) we 
compute the difference: 



/'rad,21 Prad,12 — 9 V^21 ^ 



^12 



ti=t 



ret _ j^adv 



12 



^21 



)l 



[*r'fe),t2] 



t2=t 



t2- 



Symbolically we denote 



Aret _ 
^ba — 



Aadv _ 



Vb[tf'Va)] 



(4.19) 



(4.20) 



the Lienard-Wiechert vector potential of 6-th charge at point at which a-th one is placed. 



4.3 Radiative part of mixed stress-energy tensor 

In this subsection wc find the terms which produce the radiative part of mixed momentum, 
either (4.14) or (4.15). We start with (/7-momentum (4.11) which is then nothing but the 
mixed space-time component of the following tensor: 



' int,rad 



+ 



6162 



'-hi 



_d_ 



a 



_9_ 

df2 



dt2 V ^^^dtidt2 



4 ( -Div' 



dtidtr^ 



dtidh 



(4.21) 



midt2j dt2 V dtidt2j dti 
To restore (/7-dependent terms leading to this expression, we insert 

n:,^ (/) = tr (/) - f<b (4-22) 

and substitute the right hand side of the third line of eqs.(3.16) for tt^^. (Operator TJ^'^ 
is defined by cq.(3.10).) After cancellation of like terms we obtain a linear combination 
of coefficients (3.7). Further we omit the integration over polar angle and multiply the 
result on inverse Jacobian = q/(rir2). Finally, we obtain the tensor. 



C2 



- (x,x+<xr) 



Cl 



Q3 



a 



C1C2 d'^a 
r\rl dtidt2 rj*r| dtidtl 
9V 1 



C2 



1 



a 



rfr2 dtidt2 rir2 dtfd^ 



rfr2 dtidt2 
Cl d'^a 



+ 



riri dtidtl 
1 9V 



+ 



1 



rfr| dtidt2 



(4.23) 
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which is symmetric in its indices. 

Similarly to the mixed stress-energy tensor (1.6) itself, its radiative part, 

contains the convolution of tensor (4.23), 

4 _ 9 j 2l2l . A ( _^!^^ ^ ^ A 

1 / \ C2 9V 1 9V 



r^rl dtidt2 \ dtidt2 J rirf dtidt2 rir| dtidt^ 

+ :::3:r^7^ + i2:r^;2^ • (4-25) 



rjr2 dtidt2 r\r2 dtldt2 
Integration over (p results in the combination of partial derivatives in time variables 



L 



27r 

d(fJt'^a — 6162 



a- 







d ( 


(qlkil 


dhdt2j 





dtidt2j dti \(i\\ri\\dtidt2 J dt2 \q\\r2\\ dtidt2 



(4.26) 
where expressions 

Ikall = ^JK - (van,)A;3]2 - /^2[v„n,]2, (4.27) 
a — 1,2, are introduced in Appendix A. 

Finally, after a tedious calculations^ we derive the identity d^T-^"^^^^^ = 0. It means that 
the radiative part (4.24) of mixed energy-momentum tensor density (1.6) is conserved off 
particles' world lines. 

It is worth noting that the terms in the first line of expression (4.23) only belong 
to the mixed part (1.6) of the electromagnetic field stress-energy tensor. The others 
provide vanishing of divergence of Tmt.rad as well as reasonable expression for radiated 
energy-momentum which escapes to infinity. By means of the relations 

< = la'K. K=la'{<-la'«) (4-28) 

"^a — 7a ^ai -^a ~ ^a^ai 

Ca = -ff[l + ra{ka-aa)+l-\%], (4.29) 

the expression (4.23) can be easily rewritten in manifestly covariant notations. 



4.4 Radiative part of mixed energy 



Since rf^ = —1, the convolution (4.26) contributes into zeroth component of radiative 
part of mixed energy-momentum. 



vP 

' in 



int.rad " ^1^2 

1 d 



^1^2 + 



+ 



+ 



dtidt2 
1 d 



(4.30) 



2dti \q\\ri\\ dtidt2 J 2dt2 \<i\\r2\\ dtidt2 



1 9V 



^Differentiation of Tjnt rad is straightforward: one can apply the rule (2.8) and combine like terms 
scaled as "r^"; exponents m and n run from 1 to 4 and their sum 3 < m + n < 7. 
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As could be expected the argument of mixed double derivative, 



is not smooth in neighborhoods of end points ^ E^, 



tl=t t2=t 



t2=t 
t2— >— OO 



Hence the interference rate of radiated energy heavily depends on order of differentiation, 
either d^G^ / dtidt2 or d'^G^ /dt2dti. Having used the rule (4.7) we arrive at the expression 
of type (4.14). Choosing the rule (4.8) we obtain the expression of type (4.15). The 
results can be generalized as follows: 



PUl>a= dtal-^GUtaX\ta)]+ dt aTa' G^K tf^t a)i (4.33) 

J —OO J —OO 



where 



. . J gab(^«-^&) „,-2 , QabiVg ■ '"b) , ^ 

7a <^6a = e„eb <^ -3 7(, H -3 • Vb) (4.34) 



n n 



"6 



It is convenient to rewrite eq.(4.33) in a manifestly covariant fashion, 

PUab= r dsbG:bK\sb),Sb]+ r ^^'^\sbG:b[sT{sb),Sb]. (4.35) 
J —00 J —00 

The particles' world lines Ci and C2 are parametrized by individual proper times si and 
S2, respectively. Upper limits ti and T2 label the points at which Ci and ^2 puncture the 
observation hyperplane S^. Two-point function, 

= e,e. i'^^ (1 + (fc . aJl + <^ - . (4.36) 



a a a 



is evaluated at points on the world lines of particle a and particle b which are linked 
by a null ray. Ua{sa) denotes the normalized 4- velocity of a-th particle; 4- acceleration 

< = d</dsa- 

In contrast with "one-particle" contributions to radiated energy-momentum, the mixed 
one is not uniquely defined. Having generalized the expressions (4.19) and (4.32), we 
obtain 

^•13,(1,21 ~P?ad,12 — Q2lls2^-oo ~ ^12lsi-»-oo ' (4.37) 

where 

Q2b = Qab (^ab" • ^ba)\[sa,sl-*{sa)] ' ^'^'^^^ 

This g-directed null vector is proportional to the scalar product of the Lienard-Wiechert 
potentials (4.20). 
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The difference (4.37) depends on the state of particles' motion at the end points of 
path integrals. Let us evaluate the net energy-momentum which escapes to infinity. The 
integrals over entire world lines should be substituted for the integrals over the past 
motion. We suppose that the particles are asymptotically free. Since Lienard-Wiechert 
potentials fall off at large distances inversely as the first power of the separation vector 
between the charges, the right-hand side of expression (4.37) vanishes. Hence the full 
amount of radiative energy-momentum emitted by interacting particles does not depend 
on the method of integration. 

4.5 Radiative part of mixed angular momentum 

When an one-particle problem is considered, bound components of the stress-energy ten- 
sor [3] and its torque [5] contribute into individual particle's 4-momentum and angular 
momentum, respectively. The corrections arise from the regularisation procedure which 
involves the Taylor expansion of surface integrals in which the first two terms lead to the 
diverging Coulomb self-energy and the Abraham radiation reaction 4- vector, respectively 
The finite terms depend on the form of the hole that is cut out from the integration 
hypersurface to ensure regularization. The best suited hole must be coordinate-free one. 
Teitelboim [3] integrate over hyperplane (j(t) — {ye M4 : Uf^{T){y^ — z^^{t)) — 0} which 
is orthogonal to the 4-vclocity of the charge at point z{r) at which particle's world line 
punctures o'(r). It is evident that the tilted hyperplane together with the future light 
cone cutting out the coordinate-free hole (see Refs. [3, Fig.l], [5, Fig.l], [2, Fig.5-2]). 
The considerations lead to manifestly covariant and structure-independent finite terms of 
clear physical sense. 

The difficulties associated with the computation of the mixed contribution (1.6) are 
twofold — to perform the meaningful decomposition of Tint into bound and radiative parts 
and to choose an appropriate surface of integration. The tilted hyperplane which plays 
privileged role in the one-particle radiation reaction problem is not suitable whenever 
two-body one is considered. Indeed, there is no a hyperplane which is orthogonal to 
the world lines of both the particles at all events. Kosyakov [19] constructs a piecewise 
hypersurface where a small fragment of a spacelike hyperplane E is replaced by a fragment 
of an orthogonal hyperplane (Ta{Ta) = {y e M4 : Ua^fj,{Ta){y^ — zj^{Ta)) = 0} in the vicinity 
of every intersection point. The deformed hyperplane is called locally adjusted. But the 
problem arises how to sew these fragments with E. 

Expressions (B.4), (B.5), (B.7), and (B.8) show that the surface integrals of bound 
component of the mixed part of stress-energy tensor depend on the state of particles' 
motion in vicinities of intersection points n E^. Unfortunately, the Taylor expansions 
of the divergent components of bound energy-momentum do not lead to reasonable finite 
covariant terms. It is because the integration surface E^ is tightly connected with the 
laboratory inertial frame. This choice yields the coordinate-dependent hole around a-th 
particle in the point of intersection fl E^. For this reason we assume that an intrinsic 
structure of a charged particle is beyond the limits of classical theory [2]. We do not 
require any assumptions about the particle structure, its charge distribution and its size 
(except that its "radius" does not vanish, although it is too small to be observed). To 
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reconcile the theory with observation, we assume that a dressed charged particle possesses 
finite 4-momentum and angular momentum. 

We face a problem of how the 4-momenta of interacting dressed charges depend on 
their individual characteristics such as their masses, charges, 4-velocities, etc. Valuable 
information can be extracted from the radiative part of electromagnetic field's angular 
momentum. Indeed, conserved quantities place stringent requirements on the dynamics of 
our closed system. They demand that that the change in radiated momentum and angular 
momentum should be balanced by a corresponding change in the individual momenta and 
angular momenta of dressed particles. It is felt that analysis of balance equations yields 
reasonable expressions. Indeed, it is shown [20] that the Lorentz-Dirac equation can be 
derived from the energy- momentum and angular momentum balance equations. In [21] 
the analog of the Lorentz-Dirac equation in six dimensions is obtained via analysis of 
21 conserved quantities which correspond to the Poincare symmetry of an isolated point 
particle coupled with electromagnetic field. 

Recall from Section 3 that the integral of mixed angular momentum tensor density 
over (p is the combination of partial derivatives in time variables. To reveal meaningful 
radiative terms we integrate eq.(3.26) over ti and t2 and apply Teitelboim's criteria. 

The calculation is virtually identical to that presented above, and we shall not bother 
with details. The result depends on the method of integration of mixed double derivatives, 

Mrad,a6 = r dsb{zbAG2'-ebUbAAlf) (4.39) 

J —oo 

(Symbol A denotes the wedge product.) If we choose the rule (4.8) we obtain at Mrad,2i- If 
one prefer the rule (4.7) they derive Mrad,i2- Direct calculations show that the difference 
between these expressions depends on the momentary state of particles' motion, 

M,ad,21 - Mrad,12 = ^2 A Q2l\ZZl^ ' ^1 A Ql2\ZZ^^ ■ (4-40) 

Expression (4.39) for radiated angular momentum satisfies Teitelboim criteria. In- 
deed, the integrand is finite and covariant; the radiation is accumulated with time. Non- 
covariant bound terms which are presented in Appendix B are quite different: they depend 
on the state of particles' motion at points Co H and they contain divergences. We as- 
sume that the expression (4.39) is involved in the angular momentum balance equation 
explicitly. 

4.6 Radiation of non-accelerated charges 

Let us consider a specific case of very massive particles, such that ei/mi << 1 and 
62/777,2 << 1. Electromagnetic field is too small to accelerate the charges so that the 
particles move with constant velocities. We place the coordinate origin of the Lorentz 
inertial frame at point at which the second charge is placed. If V2 = 0, the zeroth 
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component 7^ ^^21 of two-point function (4.34) is identically equal to zero. Hence the 
radiative energy, 

pL,2i= / dii7r'G'°jii,r/*(ii)]+ / dii7r'G°jii,t«*(*i)], (4.41) 

vanishes. 

If V2 = 0, the space components of this function get simplified, 

7r'G^2i = eie2^. (4.42) 

The function evaluated at advanced instant t2^^{ti) is just the function referred to the 
retarded instant t^^^iti) taken with opposite sign: G2i[^i) ^2'^''(^i)] = ~G^2i[^i; 
Because of shift in limits of the retarded and the advanced integrals the mixed part of 
radiative 3-momentum is as follows: 

PU21 = / dh^^'Gl,[h,tr{ti)]+ dti7r'G^i[ii,if^(^i)] (4.43) 

J —00 J —00 

= 6162 / dti — . 

This is vanishingly small quantity. It should be rejected if for no other reason than that 
the radiated 4-momentum can not be spacelike 4- vector. 

In a like manner we calculate the radiative angular momentum. 



^rad,21 = 6162 
^^d,21 - 0. 



r dtJtX,+"A, (4.44) 



If the upper limit of integration t — > +00, Prad2i — 0- Indeed, in this case for each 
point on Ci there are points on (^2 labeled by the retarded and the advanced instants. 
Corresponding contributions cancel each other. 

Alternative expression for mixed contribution to radiated energy-momentum. 



/+00 p+oo 
dh%'GUtT\h),h] + / dt272-'G?2[^r (^2), ^2], (4.45) 
00 J —00 



differs from 21 on the sum of null vectors. 



6162 



q[l- (nqVi)] 



ti-++cxD a 

, "'q 

t,__oo q[l + (ngVi)] 



(4.46) 

t2— >— 00 



estimated at limits of integrals. If the particles move with different velocities, they are 
asymptotically free. If both the particles are static, the difference is equal to zero because 
the distance q between charges does not change. 
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5 Equations of motion of radiating charges 



In this Section we study the energy-momentum and angular momentum balance equations. 

We calculate how much electromagnetic field energy, momentum, and angular momentum 
fiow across hypcrplane S^. We can do it at a time t + At. We demand that the change in 
these quantities be balanced by a corresponding change in those of the particles, so that 
the total energy-momentum, 

-P=Pl+P2 + ^/ dSialui{si) + ^ dS2a2^*2(s2) +Pint,rad, (5.1) 

and total angular momentum, 

M = Z1AP1 + Z2AP2 (5.2) 

2e^ /"^i 2e^ /"^^ 

+ --i / dsi {ajzi A Ml + Ml A ai) + — ^ / ds2 {alz2 AU2 + U2 A 02) + 

"J J-oo "J J-oo 



M 



int,rad ) 

"5 J-oo 

are properly conserved. 

We suppose that particles' individual 4-momenta pi and p2 are already renormal- 
ized. The words "already rcnormalized" mean that momenta contain contributions due 
to bound component of the stress-energy tensor density, including its mixed part. While 
the radiation which detaches itself from charges and leads an independent existence is 
involved exphcitly. 

In contrast with "one-particle" contributions to conserved quantities (5.1) and (5.2), 
the mixed ones, (4.35) and (4.39), are not uniquely defined. Our theory faces in the ra- 
diation problem a significant issue: non-uniqueness in determination of radiation sourccd 
by the mixed part of Maxwell energy-momentum tensor density. To solve the problem we 
mix the obtained expressions 

Pint,rad = KPrad,21 + (1 - K)Prad,12 , (5-3) 
Mint,rad = «M,ad,21 + (1 - /«)Mrad,12 , (5.4) 

and try to find out the value of constant k that accords with experience. 

Expressions Prad,a& and Mrad,a6 are based on two-point function (4.36) referred to the 
points on the world lines of particle a and particle b which are linked by a null ray, 

1 2 

<^ = -2(^1-2:2) (5.5) 
= 0. 

This implies that a displacement of 2;i(si) typically induces a simultaneous displacement 
of -22(52) because new points 2;i(si -|- 5si) and 2^2(^2 + SS2) must also be linked by a null 
geodesic, 

(^21 • Mi)dsi + {qi2 ■ U2)ds2 = 0. (5.6) 

This immediately gives 

dsa _ {qba ■ Ub) 

dSh Ta 
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where Va = -{qta ■ Ua)- 

By virtue of this expression we compare the retarded and the advanced integrals 
involved in eqs. (4.35) and (4.39). With understanding that functions Sa^"(sb) and s^^*(sa) 
are inverses, after some algebra we obtain 



ds„G^js„,.r(5a)] + g:, 



,(5.8) 



/ ds, (^, A G^f - e,«, A ) = / dSaizaAGll'-eaU^AAZ') 

J —oo J —oo 



+ 



ZaAQ 



ab 



(5.9) 



Sa— »— oo 



The relations accord with the right-hand sides of eqs. (4.37) and (4.40). 

Equipped with these relations we rewrite the mixed parts of radiated energy- momentum 
(4.35) and angular momentum (4.39) as follows: 

dSbGab[s''a\sb), Sb] + / dSaGba[Sa: sl^'^Sa)] + Qab 



Prad,ab 
M-ad,ab 



— oo 

n 



— oo 



Sa— » — OO 



+ 



d., {zb A - CbUb A Alf) 

l-Ta 

/ dSa {Za A Gll^ - eaUa A A^^) + Za A Qab 
J — oc 



, (5.10) 
(5.11) 



Sa— + — OO 



Since Qab 7^ Qba, the expressions are not symmetric in indices a and b. 

Substituting eqs. (5.3) and (5.4) into right-hand sides of eqs. (5.1) and (5.2), respec- 
tively, we obtain the total energy-momentum, 

2 



b=i I 

and angular momentum. 



2e? n 



dSbttbUb + / dSbGll + KbQba 



S;,— >— 00 



(5.12) 



^ f 2e^ r'' 

M = z]\^bApb+^ dsb {alzb Aub + UbA ab) 

6=1 1^ 

dsb {zb A Glf - ebUb A A^f ) + KbZb A Qba 

■00 



(5.13) 



s;,— »— 00 



(In our notations Ki :— 1 — k and K2 := k.) The balance equations are differential conse- 
quences of these conserved quantities. Since the action is not propagated instantaneously, 
the balance in a vicinity of the first charge as well as in a neighborhood of the second 
charge should be achieved separately. 



2el 2 

Pa — g~^a'^a ~ ~ l^aQab 

UaA [pa + -^aa - C^^^f + l^aQab 



0. 



(5.14) 
(5.15) 



27 



(The overdot means the derivation with respect to individual proper time Ta-) Solution of 
six linear equations (5.15) in four components of a-th 4-momentum contains an arbitrary 
scalar function, say nia, 

Pa = maUa- ^tta + eaAll^ - KaQab- (5.16) 

Since («„ ■ a^) = 0, the scalar product of momentum of a-th particle on its 4- acceleration 
is as follows: 

{Pa ■ aa) = + ea{Alf ■ aa) - KaiQab ' ^a)- (5-17) 

Similarly, the scalar product of particle's 4-velocity on the first order time derivative of 
particle's 4-momentum (5.14) is given by 

{Pa ■ Ua) = -fal - (G^f • Ua) - KaiQab ' ««)■ (5.18) 

Equipped with the expression (5.7) one can derive the significant relation, 

A Act 

^^a=-K-ea^: (5.19) 

where F^^ is the well-known Lorentz force of 6-th charge acted on a-th one. Substituting 
this into eq.(5.18) and summing up (5.17) and modified (5.18), we obtain 

(Pa -Ua)' ^ ea{All^ ■ Ua)' - Ka{Qab ' Ua)' (5.20) 

where dot means the a-th proper time derivative. On the other hand the scalar product 
of 4-momentum (5.16) on a-th 4-velocity is written as 

{Pa ■ Ua) = -ma + ea(A^f • Ua) - Ka(Qab ' Ua). (5.21) 

We see clearly that is of constant value. It can be interpreted as already renormahzed 
mass of a-th charged particle. 

Finally, we differentiate the expression (5.16) and substitute it for the left-hand side 
of eq.(5.14). After cancellation of like terms and taking into account eq.(5.19) we arrive 
at the Lorentz-Dirac equation, 

maaa = ^ («a - a^a) + F^f . (5.22) 

We see that particles' equations of motion do not depend on mixing parameter k. 

In terms of kinematical variables conserved quantities of our particles plus field system 
looks as follows: 

Za^ \ maUa - -^«aj ^ ~^ J A Ua + Ua A Ua) 



M 



Ta 



J2 I dSaZa A F^^\ (5.24) 
a^b 
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The work done by Lorentz forces of charges acting on one another exhausts the radiation 
reaction due to combination of fields. 

Individual 4- momentum (5.16) of a-th dressed charged particle is modified comparing 
with the well-known Teitelboim's expression (1.7). The bound component ilnt.bnd of 
mixed part of the electromagnetic field stress-energy tensor contributes two additional 
terms: "immovable core", CaAl'^^, of clear physical sense and "changeable shell", —KaQab, 
that heavily depends on mixing parameter k. The corrections are inspired by unavoidable 
deformation of bound electromagnetic "clouds" due to mutual interaction between the 
sources. In my opinion, the changeable term arises due to forced choice of non-covariant 
surface of integration (see Subsection 4.4). For this reason we proclaim the expression 

Pa = rriaUa - -ytta + CaA^^ (5.25) 

as the only one of true physical meaning. 



6 Discrete symmetries 

In the past Sections we have emphasized the importance of the invariance of the action 
integral (1.1) under the continuous group of space-time translations and rotations. Ac- 
cording to Noether theorem, these symmetry properties imply conservation laws, i. e., 
those quantities that do not change with time. In this Section wc study symmetry prop- 
erties of energy-momentum and angular momentum carried by electromagnetic field which 
rely on invariance of (1.1) under discrete transformation groups. 

6.1 Time inversion 

The transformation of the time inversion T is defined by [2] 

T:y''^y'^ = y^. (6.1) 

It immediately gives 

T:^,^(sJ^y^«)=^„,^(s„). (6.2) 
The proper time possesses odd parity [2], 

T : s„ 1-^ s'„ = -s„. (6.3) 

The other kinematic quantities then follow easily, 

r : u>^^{sa)^u':{s':)^-UaASa), (6.4) 

r : <(s„)^at«)=a„,^(s„), 
etc. The retarded and the advanced instants transform into each other, 

rC{s,) = sf^s.l Tsl'^is,) = sT{s,). (6.5) 
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Inserting eqs. (6.4) and (6.5) into two-point function (4.36) yields 

rGX\s,),s,] = ^Ga,A^f^{s,),s,], (6.6) 
TG^Js^ (.,),.,] = -GaUs:^\s,),s,]. 

To establish the symmetry properties of Prad a-nd Mrad with respect to time inversion, 
we locate the observation hyperplane at the distant future, 

^2 /• + CXD 0^2 /• + C30 



rad 



rad 



2g2 /•+00 2g-^ f+°^ 

= — i / dsia^Mi + ^ ds2alu2 (6.7) 

/+00 i>+oo 
ds, {Glf + G^f ) + (!-«:)/ ds2 (G^^' + ) , 
oo J —oo 

= --i / dsi (a^2;i A Ml + Ml A ai) + — ^ / ds2 (^2^2 A M2 + M2 A 02) 
3 J-00 3 J -00 

/+00 
dsi [^1 A + Gt^-) - eiMi A {Af^ + )] (6.8) 
-00 

/+00 
dS2 [Z2 A (Glf + G^^) - 62^2 A (^llf + ^if )] ■ 
-00 



These expressions give the full amount of radiation emitted by interacting particles. 

The transformation of the retarded and the advanced functions Gab into each other 
implies that the radiated energy-momentum and angular momentum are of odd time 
parity, 

'^^ad = -^rad,M , ^M^ = "^rad,,. ■ (6-9) 

6.2 Space inversion 

This operation is defined by [2] 

V:y^^y'^^ -y^. (6.10) 

Correspondingly, 

P:<(.„)^z'^«) = -z„>„). (6.11) 

The proper time Sa remains invariant. These transformation properties imply that that 
particles' 4-velocities and 4-accelerations change as follows: 

r : <(s„)^<«) = -M„,^(s„), (6.12) 
V : a^(s„)^a'^(s;) = -a„,^(s„). 

Since the retarded and the advanced instants remain invariant, the basic two-point func- 
tion (4.36) transforms analogously, 

= -G„,,„. (6.13) 
Substituting this into eqs. (6.7) and (6.9) and using the relations (6.11) and (6.12) yields 

VP,^^ = -P.ad,M , VK:^ = Mrad,^. . (6.14) 
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6.3 Reciprocity of particles 1 and 2 



When a closed system of two identical charges is considered, radiative conserved quantities 
should be symmetric in indices 1 and 2 that label the particles. Having interchanged these 
indices in function (4.36) we obtain 

G^dsi,sr{s^)]l^^ = GUs?\s2),S2], (6.15) 

From the reciprocity relations we see clearly that mixed parameter k in expressions (6.7) 
and (6.9) should be equal to 1/2. Choosing the linear superposition 

Pint.rad = ^ (Prad,21 + Prad,12) , (6.16) 
Mint,rad = ^ (^M,^d,2l + M,^,^), (6.17) 

we restore invariance of radiated energy-momentum and angular momentum with respect 
to reciprocity of particles 1 and 2. 



7 Conclusions 

The present paper is devoted to study of phenomena of emission and propagation of 
energy in classical electrodynamics. The field in action (1.1) has its own uncountably 
infinite degrees of freedom. Variation of (1.1) yields Maxwell's equations with point-like 
sources and equations of motions of particles interacting through the medium of the field. 
The problem then becomes one of mutual determination: the field is determined by the 
charged particles and their motion, and the motion of the charges is determined by the 
field. 

In this paper we study interference of outgoing electromagnetic waves in a hyperplane 
T,t = {y E M.4 : = t} associated with an unmoving inertial observer. We calculate how 
much electromagnetic field energy, momentum, and angular momentum flow across this 
hyperplane. Surface integration of the stress-energy tensor (1.4) over reduces field's 
uncountably infinite degrees of freedom. After the renormalization procedure we arrive at 
the action at a distance theory [22, 23] where particles interact directly with one another. 
The fields in resulting expressions (5.23) and (5.24) do not have degrees of freedom of 
their own: they are functionals of particle paths. Following Ref . [24] , we refer to them as 
direct particle fields. 

Starting with the retarded Lienard-Wiechert fields, after integration we arrive at the 
retarded and the advanced direct particle fields. The retarded and the advanced instants 
arise naturally as the end points of interference integrals (2.15) and (2.16). But the 
retarded causality is not violated because the advanced instants as well as the retarded 
ones are before the fixed observation moment t. Direct particle fields referred to advanced 
instants do not describe neither incoming radiation nor converging electromagnetic waves. 

Nevertheless, the retarded and the advanced quantities transform into each other under 
the infiuence of inversion of time axes. To demonstrate the invariance of electromagnetic 
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field's energy-momentum and angular momentum with respect to time inversion, we locate 
the observation surface at distant future. We show that the full amount of radiation 
emitted by a closed system of two interacting charged particles is invariant with respect 
to past and future as well as with respect to inversion of space axes. 

The part of Noether quantities which escapes to infinity is defined by basic two-point 
function (4.36). It is equal to the sum of Lorentz force Fab of c^-th charge acted on 6-th 
one and the total time derivative of direct Lienard-Wiechert potential Aab, taken with 
opposite sign (see eq.(5.19)). In the specific case of very massive particles, such that 
ei/mi « 1 and 62/777.2 << 1, the mixed energy- momentum Pint,rad vanishes as could be 
expected for nonaccelerated charges. Indeed, let us consider static charge 6 at a coordinate 
origin. Zcroth component F^"^, either retarded or advanced, is potential one: its integral 
over indicated portion of (a cancels the change of only nontrivial A^^. Space components 
F^g^ of the retarded and the advanced Lorentz forces compensate each other. Having 
performed a trivial Lorentz transformation we extend the statement on a motion with 
constant velocity. Radiative angular momentum possesses analogous properties. 

To derive the equations of motion of interacting particles we compare fiows of energy- 
momentum and angular momentum through very close hyperplanes and ^t+At- Having 
balanced particles' individual characteristics and corresponding quantities carried by elec- 
tromagnetic field, we obtain the well-known Lorentz-Dirac equation of motion of charged 
particle in the retarded field of the other charge where the self-action is taken into ac- 
count. Since the Lorcntz-Dirac equation possesses pathological solutions, such as runaway 
solutions or preaccelerations, the authors [25, 26, 27] propose the Landau-Lifshitz equa- 
tion [9, §76] as more satisfactory alternative. Spohn [25] showed that a solution of the 
Lorentz-Dirac equation which does not satisfy the Landau-Lifshitz equation is of the run- 
away type. Moreover, Landau-Lifshitz equation does not permit runaway solutions or 
preacceleration [27]. 

Obtained expressions for radiated energy-momentum and angular momentum (includ- 
ing interference effects) are valuable also for a strictly prescribed motion of particles under 
the influence of a very powerful external force. It is necessary to compare them with cor- 
responding results for circhng charges [12, 13]. It would be interesting to consider the 
specific case when the charged particles are a very close to each other. Since the electro- 
magnetic field satisfies the superposition principle, the models either an extended object 
consisting of point charges or a continuous charge distribution are based on dynamics 
of two-body system [28, 29]. 
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Appendix A Integration over angular variable 

To calculate the mixed rates of energy-momentum (2.1) and angular momentum (2.2) 
carried by the electromagnetic field, we should first perform the integration over angle. 
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When facing this problem it is convenient to mark out 99-dependent terms in expressions 
under the integral sign. In the Maxwell energy-momentum tensor density we distinguish 
the second-order differential operator (3.10) with (/^-dependent coefficients (3.7). It can 
be decomposed into a combination of partial derivatives in time variables 11^ given by 
eq.(3.13) and tail tt^ of the type in eq.(3.14). 

This Appendix is concerned with the computation of the tails. Equipped with them 
we express the aforementioned integrand as a combination of partial derivatives in ti and 

t2- 

To implement this strategy we must first integrate the coefficients (3.7) over the angle 
variable. We start with the simplest one, 



f2TT 

qrir2 



where numerator a is equal to 1 or K^. Our task is to rewrite the integrand as a sum of 
term with denominator ri and term with denominator r2- To do it we introduce a new 
layer of mathematical formalism and develop convenient technique. 

Let V be the vector space such that io, ii, and '12 is its hnear basis. We shall use rfap — 
diag(— 1, 1, 1) and its inverse ■q"'^ — diag(— 1, 1, 1) to lower and raise indices, respectively. 
We introduce the pairing 

(■) : V xV (A.2) 
(a ■ b) ^ r?„^a"b^ 

which will be called the "scalar product" . 

Wc introduce null vector n" = (1, sin 99, cose/?) which belongs to the vector space V. 
We express the (^-dependent constructions 

ra^-{Ka-Va), Ca ^ + (Ka ■ Va) (A.3) 

as the scalar products —{r^-n) and (Ca • n), respectively. We shall use sans-serif letters for 
the components of timelike three- vectors ra eV and Ca e V, 

r° = - (v„nq)A:f, ^l = hvaW, = /^^t'aW; (A.4) 

Ca = -7a"' - (Vanq)A;^, ^ hVaW , 4 = hVaW , (A.5) 

where orthogonal matrix, 

(cos(/9^ — sin(/?g \ / cos-i^g sin-i^^ \ 
sin^jg cos (pq 1 , (A.6) 

1 / \ -sint?5 cos^?5 / 

is constructed from components of the relative position 3- vector q = Zi — Z2. The numer- 
ator in eq. (A. 1) is the scalar product (a • n) = — a° + sin ip + a'^ cos (p; it is equal to 1 if 
vector a = (—1, 0, 0). 
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We introduce the dual space of one-forms, say W, with basis lj^,u!^ and a)^, such that 
= 5^, where io, ii, '12 constitute the basis of V. The wedge product L = a A 6 of two 
one forms a and b constitutes two-form, 

L = {aobi - aibo) A (2;^ + (ao&2 - 0-2^0) 00° Aui^ + (0162 - 02^1) A u^. (A.7) 
We introduce dual three-vector L —*L with components 

r = ^e-^'L^, (A.8) 
= e^'^^apb-y. 

^oiPi denotes the Ricci symbol in three dimensions, 

{1 when q;/57 is an even permutation of 0,1,2 

— 1 when a(5^ is an odd permutation of 0,1,2 (A. 9) 

otherwise. 

We raise indices in eq.(A.8) and define the vector product of two vectors, a and b, 

r = sVa'^b". (A.IO) 

Tensor, 

= e'^^^V^t^V^u. (A.ll) 

has the components 

Now we calculate the double vector product, 

[A[BC]]" = sVA^^VB'^C^ (A.13) 

Since 

= -<^%?7/3. + (^".r?^^, (A.14) 

we arrive to the unusual rule, 

[A[BC]] = -B(A- C) + C(A- B), (A.15) 

instead of the well-known law acting in space with Euclidean metric. 

To simplify the denominator rir2 in the integrand of cq.(A.l) as much as possible, we 
rewrite 27r-pcriodic functions Ta = —^afi ~ ^a,i sin 99 — ra^2 cos 99 as follows: 



ra = -ra,o - Pa sm{(p + Pa = ^ r^i + r^,2- (A. 16) 

(We recall that is the scalar products (r^ ■ taken with opposite sign, components 
are given by eqs.(A.4).) Shift in argument of harmonic function is determined by the 
relations 

fa,! = Pa COS 0„, r„,2 = Pa sin 0„. (A. 17) 
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After some algebra we rewrite the integrand of eq.(A.l) as the following sum: 
a + cos((^ + 0i) Al-^ - C1^p2 cos{ip + ^a) 



+ 



rir2 ri r2 

where a — (a • n). Coefficients 74^2, ^21 ^12 satisfy the vector equation 

—-^12^2 — -^21^1 + C'i2Li2 = a, 

where by ri and r2 we mean three- vectors with components in eq.(A.4) and L12 = 

To solve equation (A. 19) we postmultiply it on the vector product [riLi2], then 
vector product [r2L2i], and, finally, on L12. After some algebra we obtain 



(A.18) 



(A.19) 
on the 



^12 



([ari] • L12) _ ([ara] • L21) _ (a • L12) 

J^12 J^21 J^12 

where the denominator D12 — {\-i2 • ^12) is symmetric in its indices. 
Substituting eq.(A.18) into eq.(A.l) and using the identities 



1 

2^ 



27r 



dip 



Pa Sm{ip + (pa) 



1 



1 



Pi 



(A.20) 



(A.21) 



1 ^ COSI 

Jo 'a Pa 



COs{(fi + (f)a) 



sm{(p + (f)a) 



yields 



pa 



Aa Aa 
^12 _|_ ^21 



q||ri|| q||r2 
after integration over ip. 

Now we turn to the calculation of the coefficient 

-.271 

j^a 



27r./n qri(r2)^ 



(A.22) 



(A.23) 



Equipped with the relations in eq.(A.18) wc rewrite the integrand as a sum of terms which 
are proportional to the 1/ri, l/r2, and l/(r2)^, respectively. Using the identities 



1 r^"" 

^Jo K - Pa 



2-K 



sm{(p + (f)a)f 
C0S((/? + (f)a) 

SttJo "''[rO-p„sin((^ + 
and taking into account the relations in eq.(A.21) gives 
1 



(A.24) 



= 



+ 



W2 
1 

qlkill 



^^12^12 



D21 ' ■ q||r2| 
([ari] • [C2ri])" 



AC2 Aa I Aa AC2 _ 
^12^21 ^ ^19^91 



(a ■ C2)(ri • ra) 



12^21 



21 



D 



12 



(A.25) 
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The resulting expression for the term 

pa ^ 

^ "27r 



2tv 



dip 



aci 



q(ri)2r2 



(A.26) 



can be obtained by interchanging indices 1 and 2 in the right-hand side of eq.(A.25). 

After a routine computation based on the repeated usage of relation (A. 18), we find 
the most complicate term, 

"2^ 



1 



d(p 



aciC2 



q(ri)2(r2)2 



(A.27) 



where 

7 _ 9401402 ([ciri] ■ [C2ri]) 



B12 B21 (a • ri) (a • ra) 

"I n ^7 + >^l~ii iTq + 'J2~ ^75"' 



qiK2i 



qiKii 



qiK2i 



D 



(A.28) 



12 



12 



"^^12^12^21 "J^l 9^1 9^91 -^^19^19^91 \ 



■^12^12^21 



12^12^21 



■■12 



([ciri] • [csrs]) , ([cira] ■ [csri]) 



^21 



-D12 -(^12 

ci ([C2ri] ■ [ari]) ([ciri] • [ari]) , ([cafi] ■ [ara]) , ([ciri] ■ [ars]) 

n 21 r) 12 +^12 , 

-t>'l2 -t^l2 -t^l2 -t^l2 



+ 



D 



and the others, B21 and J2, can be obtained via interchanging indices 1 and 2. 

Wc now turn to the differentiation of coefficient (A. 22) with respect to time variables 
ti and ^2- Having substituted o = (—1,0,0) for a in the expressions (A. 22), (A. 25), and 
(A.27) we obtain the terms T)^. and J\P, respectively. The remaining term, C*^, can be 
obtained from via reciprocity. The calculations are based on the relations obtained 
via differentiation of third components and h?, i.e. the square of radius of the circle 

c{o,h)^s^nS2, 



dta 

dti 

They immediately give 



-1)"^ + (n,v„) 



q 



-2k 



q 



ml 



2ki 



3 ' 2 



q 



dta 



-11 



q 



:v„nq)r° + [v„nq]2A;f 



(A.29) 



(A.30) 



(v„nq)r° + ([v„nq][v6nq])A;^ 



and eventually give 
d{ra ■ fa) 



dta 

djra ■ ra) 

dt, 
djrb ■ rg) 

dt. 



= 2(-l) 



(ra-0 + (-l)«^^(r„.r„) 

K rg(ra ■ r„) - rO(r„ ■ rfe) 



(A.31) 



q 



- {ra-c,) + {-lf 



q 



q 



Jb ■ ra) + 



(vaHq) _ ^ kl rO(r6 • ra) - rO(r6 • r^) 



q 



q 
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We use Latin indices a and b which run from 1 to 2 (a 7^ 6). We use bold script for 
conventional three- velocities = dzl/dta and unit three- vector riq = q^q in q-direction. 
By (vfeHq) and [v^nq] we denote the conventional scalar product and cross product of 
these vectors, respectively. 

Usage of these relations allows us to calculate the derivatives of coefficients 

^ = + + (A.32) 

Ota J^ab 

^^ab _ r^A^bAO ([org] • [CferJ) I r°.1° -^r°/1°'l 



q q 
Substituting these into equality 



dta dta Vqlkill q|l''2 

and using the identities 



«'°=«f4L + 4Ll (A.33) 







dta \ 








dth \ 





yields 



= _ 1 Klo) , (A.34) 

\Kr Fall q 

\ ^ / ^b (ra ■ rfc) (nqVg) fcg r°(r„ • r^) - r°(r„ • r„) 



q^ll''2||^ ' dt2 q^lkill^ 



Further we calculate the partial derivative dC^/dt2, subtract it from A^, and prove the 
identity 

(One can derive dB^/dti, subtract it from A^, and compare the result with 9/9^2 (C° — 

dvydti).) 

Now, we calculate the tail 



where 
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The zeroth component, = k^, does not depend on ip. Inserting relations 7?° = k'^V^, 
B° = k'^B^, Cl = /c°C°, and ^° = k'^A^ into eq.(A.37) and taking into account identity 
(A.36) yields 

Space components, K^, depend on ip. They can be expressed as the scalar product 
(K^ • n) where components of three- vectors K\^V are as follows: 

Kfi^<kl Ki,i = Wi, K,2-hw'2. (A.40) 

Here uj'^j are components of the orthogonal matrix (A. 6). Having substituted KJ^ for a in 
expressions (A. 22), (A. 25), and (A. 27) we obtain the terms P*, and A\., respectively. 
The last term, Q, can be obtained from B\ via reciprocity. To differentiate them we need 
the equalities 

= (K..cO-.;r?-<^^(K;.rO-^(r..r.) 

011 q q 

- ^ [Kl,o(r, • rO + r?(Kl • n)] , 

- -^ir^ - (K, • rO - -(ri ■ r^) 

+ ^ [^2(K1 • ri) - Kl,o(ri • r^) - 2r?(Kl • r^)] , 
— ^7 ^ = (Kl • r2 + ^ ri • r2 

q q 
+ ^ • '•2) + K;,„(r, . „) + 2r°(Ki • r,)] , 

= (K;.c.) + (^^a^(K;.r.) + !^(r..r.) 

012 q q 

+ ^ [Kl,o(r2 • r^) + r°(Kl • r^)] , 



in addition to eqs.(A.31) and (A. 34). 
The derivation of equalities 



- (Kl,o[n,vi]^ + r?[n,[vin,]]' 
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is virtually identical to that presented above, and we shall not bother with the details. 
By virtue of the relation K2 — K\ — the integrals (A. 38) are related as follows: 

Vi = V\ + q'V^, B\ = B\ + q'B^, C^ = Cl + q'C^, Ai = A\ + q'A^. (A.42) 

Equipped with these expressions we find 

Q-fi . 1 Kza)(K..,) (A.43) 

oti q II ''2 II q||''2||'^ q 

dt2 II rill q||ri||'^ q 



q||ri||3 q 



Finally, after a straightforward (but fairly lengthy) calculations we derive the following 
relations: 



which generalize eqs.(A.39). 
We will need also the tail 



where 



a/3 _ O -^12 0'Pl2 0'^12 , .a/3 ( h 



-n"/3 /"^"j ^1^2 oa/3 1 Z^'' , K^K2C2 . 

,a/3 _ f'^ , K-KIc, ,a/3^^ ['^,.^^11^1^ 

\2- 



27rJo q(ri)V2 27r Jo q(n)n^2) 

It can be obtained by means of covariant generalization of previous relations. Setting 
a = and (5 — Q and taking into account eq.(A.36), we obtain 



TT 



00 _ 7.0 / /70 
12 



^0 / ^0 _ +^0 /'^O _ ^ j ^pO 



dti J '\ dt2 

where relations "D^ = k^T^^-, ^1 = kl^B^, and C° = A;°C° are taken into account. If a = i 
and /3 = 0, we have 



^ + .^ (b^-^ 



- Cl-^ + vlfB^.-^^-vlV'^. (A.48) 
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If a = and P = j, we arrive at 



dto 



An obvious generalization of expressions (A.47)-(A.49) is 



TT 



al3 
12 



(A.49) 



(A.50) 



Since the angular integration leads to a combination of partial derivatives in time 
variables, the end points are valuable only. At these points the radius of circle C(0, K) = 
Si n S2 vanishes. Hence we can restrict oneselves to calculation of expressions — c?2^^i2 
and — diVi2 at points where h = 0. 

To simplify the calculations as much as possible we express the integrands of eqs.(A.46) 
in form of expansions in powers oi h. It allows us to remove harmonic functions from de- 
nominators. Since the derivatives dh'^/dta do not vanish whenever h!^ — (see eq.(A.29)), 
it is sufficient to expand the first order of this parameter, 



^12 



1 + 



nnqVi]2 ^ [nqV2]' ^ ([nqVl][nqV2]) 



M\2 



M\2 



.0^0 



1' 2 



(A.51) 



+ 



+ 



+ 



/i2 kfn\[n^[yin,^]y + A;fn^[nq[vinq]]^ 



0^2^0 



2 q(rO 
/cf<[nq[v2nq]]^' + k^n{[nq[Y2nq\Y 



qr?(r0)2 



J^2 ^tj _ 

2 qr°r2 



With degree of accuracy sufficient for our purposes. 



B 



12 



C 



Kkln{klcifl 



(A.52) 



q(r°)^r2 

The calculation is straightforward, although it involves a fair amount of algebra. Finally, 
we obtain 



*-'l2 



^^3 
'-'12 



dto 



dU 



kfnl^ 



^q^l^q^2 



kX { CI 



n^kinljv2 I C 



+ khv' 



du 




q2j,o 



(A.53) 
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Appendix B Bound parts of mixed energy- momentum and angular momen- 
tum 

In this Appendix we consider tlie parts of mixed energy and mixed momentum which 
describe unavoidable deformations of electromagnetic "clouds" of charged particles due 
to mutual interaction. The short-range terms will be absorbed by four-momenta of "bare" 
particles within the renormalization procedure. They arise from the bound part of mixed 
momentum, 




which is equal to the difference of the total interference momentum (3.18) and the radiative 
component (4.11). 

Our next task is to integrate Pfnt.bnd ^^^^ time variables by means of expressions (2.15) 
and (2.16). Recall from Section 4 that it is sufficient to evaluate the arguments of time 
differential operators at the ends of integration intervals. 




Figure 6: Boundary points of interference of spherical wave fronts 5*1 and 5*2 in a hy- 
perplane S^. In momentarily comoving Lorentz frame pictured in Fig. 1 the wave fronts 
intersect at the coordinate origin, li h — 0, triangle OiO^H reduces to the hne. Distance 
q between particles is equal to the difference of radii of spheres. 

In Figs. 6 the boundary conditions of interference of spherical wave fronts 5*1(01, A;^') 
and 52(02, /C2) are presented. (Fig. 3 pictures the combination of waves in four- dimensional 
spacetime.) Distance q between their centers, Oi and O2, is equal to the difference of their 
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Figure 7: Acausal interference in a hyperplane E^. Distance q between particles placed 
at centers of spheres is equal to the sum of their radii. 



radii, and A;^. In the left figure q = k2 — while in the right one q = A;° — /c^- In- 
serting this into eq.(2.5) gives = A;° and k^ = —k^ for the first and the second cases, 
respectively. 

The retarded and the advanced instants label the points on particles' world lines which 
are connected by a null ray. Since the relative position vector q is of null length, the null 
vectors Ki and K2 are coUinear: (Ki ■ K2) — 0. Their space parts Ki and K2 are 
codirectional in this case. 

Within acausal region wave fronts combine in quite different manner (see Fig. 7). In 
contrast to the pair of points Ziiti) and Z2[f2^{ti)], the vertices Ziiti) e Ci ^'^^ ^2(^2) ^ C2 
are spacelike related. Indeed, the scalar product of the separation null vectors takes 
minimal value —2k\k2 if their space parts Ki and K2 are opposite directed. Hence 



In Table 1 we collect the basic quantities and functions which are evaluated at the 
boundary points. 

To integrate the expression (B.l) over time variables, we apply the scheme developed 
in Section 4. Usage of the rule (4.7) implies 




6162 V2 - n; 



Pbnd,21 — 



+ 



2 A;0 [1 + (n,V2)] 



(B.2) 



+ 
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Table 1: Relative position, retarded distances and some basic functions at the ends of 
integration intervals 





[tT\t2),h], [ti,tr(ii)] 


[h,t'^{t,h)], [t[{t,t2),t2] 


q = k^- kl, g° = +q 

1,3 — UO 1.3 _ uO 

rv-^ — "^l 5 '^2 — 2 

r^ = A;? [1 - (n,vi)] 
r2 = A;0 [1 - (n,V2)] 


q = A;'j' — A;2, g'^ = — q 

l3 _ _i,0 p _ _J,0 

rv-^ — ''-'1 5 ''^2 — 2 

n = A;? [1 + (n,vi)] 
r2 = A;° [1 + (n,V2)] 


q = A;? + A;0 

P — —h'^ — P 

rt-|^ — rv-^ 5 rt'2 — "'2 

n = A;? [1 + (n,vi)] 
r2 = A;0 [1 - (n,V2)] 


a = 
^ = q[l-(n,vO] 

|| = -q[l-Kv2)] 


CT = 
|| = -q[l + (n,vO] 

^ = q[l + (n,V2)] 


a = -2A;0A;0 
|^ = -q[l + (n,vi)]+2fc0 

|| = -q[l-Kv2)] + 2A;? 


Ai = A;?q [1 - (n,vi)] 
A2 = -A;°q[l-(n,V2)] 


Ai = -A;Oq[l + (n,vi)] 
As = A;Oq [1 + (ngV2)] 


Ai = -A;S'q[l + (n,Vi)] 
A2 = -A;0q[l-(n,V2)] 



If one prefer another order of differentiation (4.8) they obtain 



Pbnd,12 



6162 v{ + nl 



+ 



2 A;? [1 - (n,vi)] 
6162 v\ + ni 



+ 



6162 v{ - n; 



2 A;?[l + Kvi)] 



[q-*{t),t] 



t2— + — 00 



(B.3) 



The lower limits ta ^ —00 vanish even if the motion is finite. Final expressions depend 
on particles' positions and velocities referred to the moments and t'^2^{t) as well as 

on the laboratory time t itself, 



Pbnd,21 



6162 



[1 - (n,V2)] 



[t>*r*(*)l 



_|_ 



2 t2-t A;^ [1 + (ngV2)] 



6162 ,. 
hm 



Vn — n„ 



2 t2-t A;0 [1 - (ngV2)] 



(B.4) 



Pbnd,12 



6162 



6162 



v\ + n* 



2 A;? [1 - (n,vi)] 



[ti,t-='(ti)l 



lim 



2 A;? [1 + (rigVi)] 



[tl,t'2(t,t2)] 



(B.5) 



In an analogous way we find short-range contribution to the mixed energy due to time 
integration of the following expression: 



-po — T>o _ 7?o 

' intjbnd ' int ' int,rad 



(B.6) 



= £162 
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1 d 



9V 



1 d 



2 9^1 VQlkill ^^i^^2/ 2 9^2 \q| I ''2 1 1 

The calculation is virtually identical to that presented above, and we shall not bother 
with details. Finally, we obtain 



Pbnd,21 



Pbnd,12 



[1 - (n,V2)] 

6162 



k1 [1 + (n,vi)] 



[tT'{t),t] 



+ hm — Q- 

t2^t k2 



1 



1 + (ngV2) 
1 



- 1 



[«r*(*2),t2] 



(B.7) 
fB.8) 



As could be expected for bound terms, they (i) depend on the momentary state of particles' 
motion, (ii) contain divergent terms, and (iii) are non-covariant. 

The bound components of angular momentum have similar structure. 



-"^^bnd,21 



+ 



6162 

6162 



t4 + 



^2% 



nM5=*W] 



'^2°[1-KV2)] fcO [1 - KV2)] 



(B.9) 



lim 



v\ + n' 



^^2°[1 + KV2)] kl fc2°[l + (n,V2 



t- 



n„ 



2 t2^t 

6162 

lim V 

2 t2^t \_kl[l- (riqVa)] k^ 



^2 



^'^-'bnd,21 



+ 



6162 
6162 

2 

6162 



1^2 



A;" [1 - (n,V2)] 



(B.IO) 



+ 



lim 

2 t2-^t 



^2l^^2 



1 + (ngV2) 1 - (nqV2) 

[<r'(t2),t2] 



n„ 



fco [1 + (n,V2)] 



6162 ~^v-:^ 
hm 

2 t2^t 



[1 - (n,V2) 



[t[it,t2),t2] 



Alternative expressions, Mj^^^^^, can be obtained via reciprocity of indices 1 and 2. 

In contrast to one-particle case, expanding of the expressions under limit signs in 
powers oi Aa — t — ta does not yield simple and manifestly covariant terms of clear 
physical sense. The "deformation" is due to the choice of the coordinate-dependent hole 
around the particle in the integration surface E^. We neglect these structureless terms. 
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